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ABSTRACT 

Stars form with a complex and highly structured distribution. For a smooth star cluster 
to form from these initial conditions, the star cluster must erase this substructure. 
We study how substructure is removed using N-body simulations that realistically 
handle two-body relaxation. In contrast to previous studies, we find that hierarchical 
cluster formation occurs chiefly as a result of scattering of stars out of clumps, and 
not through clump merging. Two-body relaxation, in particular within the body of 
a clump, can significantly increase the rate at which substructure is erased beyond 
that of clump-merging alone. Hence the relaxation time of individual clumps is a 
key parameter controlling the rate at which smooth, spherical star clusters can form. 
The initial virial ratio of the clumps is an additional key parameter controlling the 
formation rate of a cluster. Reducing the initial virial ratio causes a star cluster to 
lose its substructure more rapidly. 
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1 INTRODUCTION 



The vast majority of stars do not form alone. 
They appear to form in a hierarchy in structures 
of t ens to tens of t ho usand s of stars jTesti et al.l 
2000l: iGutermuth et al] l2005l: ISanchez et all l20o'7f 



In a clustered 
may occur such as 



Andre et al.l l2007l: [Goldsmith et all l2008l: lAndre et all 



20101; iBressert et all |2010|: Idi Francesco et al.l l2010l : ^Hn^m^^^ 
Gutermuth et all l2009l : iGouhermis et all I2OI0I ). Such ° ^ a 



phase many interesting processes 

„ rap id dynamical mass segregation 

jAllison et al.l 20091. 2010l). binary disruption and modifica- 
tion (|Heggidll974l : lKroupalll995l : [p"arker et al.ll2009t ). the for- 
mation of high-order multiples like t h e Trapezium system 
jAarseth et al.l 1 1974 IZinneckej I2OO8I : lAllison et all |201G| '). 
and s tar-disc interactions affecting planetary system forma- 

Pfalzner et al] 



structures are a natural conseque nce of the gravo-turbulent 
model of star formation (e.g. [ Klessen fc BurkertI I2OOOI: 
Bonnell et al.l l200ll. lioosl: iBate et al.l l2003l : iBonnell et al l 
20081 : lBatell2009l : IOffner et al.ll2009l ) 

Hierarchical distributions are not in equilibrium, and 
will rapidly dynamica lly evolve into dense star clusters or 
loose associations (e.g. Aarseth fc Hillslll972l: lGoodwin!ll998" 
Bate et al.|[l998l; iBoilv et al.lll999l: iKroupa fc Bouvicr 20^ 



I2OO5I : iThies et al 



1998 



2005 



Watkins et al][l99; 
2OIOI I 



Therefore an understanding 
of the collapse of hierarchical distributions is important to 
understand the formation of star clusters and the possible 
importance of these effects. 

The evolution of initially substructured stellar 
distributions into smooth star clusters has been stud- 



Good win fc Whitworth' '20 041: lAUison et al.' 200 9.. I201ol: 



Moeckcl fc BonncU 2009 : Fellhauer et al.| [2009: Gi cles et all 
2OIOI ). Such evolution is espe cially violent if the stars are 



ied by many authors (e.g 


. Aarseth fc Hillsl 19721: 


GoodwinI Il998l: iBoilv et al. 




I1999I: Kroupa fc Bouvierl 


20031: Goodwin fc Whitworth 




20041: Allison et al. 12009: 


Fellhauer et al. 20091: Moeckel fc Bonnell 20091: Smith et al.l 


201lf). In particular, iFellhauer et al.l (|2009f) attempted to 



initially dynamically cool (see lAUison et al.|[2009l. |2010|) as 
many observations^ sugge s t they are (e.^ . Walsh et al. 20041: 



Di Francesco ct al. 2004; Pcrctto ct al. 2006: Walsh ct al 
2007 ; Andre ct al..,2007 ; .Kirk ct al...2007. : .Gutermuth ct al 
200i ■ 



E-mail:rsmith@astro-udec.cl 



quantify how rapidly an initially clumpy distribution could 
evolve into a smooth star cluster. 

In this paper w e p articularly extend the work of 
IFellhauer et all (|2009l ) and lSmith et al.l (120111 ) to investigate 
how a collapsing clumpy distribution in a static gas poten- 
tial is able to erase i t s sub structure and form a smooth clus- 
ter. IFellhauer et al.l (|2009l ) presented a semi-analytic model 
for the erasure of substructure, but they did not properly 
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account for two-body effects in their simulations. Here we 
revisit their analysis with an accurate A'^-body code. 

In Section 2 we present our initial conditions, in Sec- 
tion 3 we study the erasure of substructure before examin- 
ing the rates at which substructure is erased in Section 4. 
Finally we discuss our results in Section 5, and draw our 
conclusions in Section 6. 



2 INITIAL CONDITIONS 

We perform our A^-body simul ations using t he direct A'^- 
body integration code nbody6 (|Aarsethll2003l 'l. The advan- 
tage of NBODY6 is that it is able to rapidly and accurately 
model stellar dynamics, and two-body encounters in partic- 
ular. 

Our initial co nditions are similar to those of 
iFellhauer eraLll2009l (hereafter F09). Our young star form- 
ing regions have a total mass of lOOOM©. We assume they 
convert gas to stars with an efficiency e which ranges be- 
tween 0. 1 and 0.8. Therefore the mass of stars is e x IOOOM0 , 
and the mass of gas (1 — e) x IOOOMq. 

We simulate the gas with a static background Plummer 
potential with a Plummer scale radius of R^i which ranges 
between 0.02 and 1 pc. We set a limiting cutoff radius for 
the gas potential of 5 x flp"}. We acknowledge that this is 
not ideal, as the gas will also dynamically evolve. However, 
it is presently impossible to simply model live background 
gas and so we f ollow previous studies in inclu ding a static 
background (e.g. lMoeckel fc Bonnelll2009l , F09 ISmith et all 
[20il 

The stars are distributed within the gas potential in A'o 
subclumps which follow the underlying gas Plummer dis- 
tribution. A'^i) ranges from 4 to 32 resulting in a mass per 
clump of Afpi = 6 to 8OA/0 (where the clump stellar mass 
is A/pi = (e x lOOO)/iVoM0). 

Subclumps are distributed within the Plummer sphere 
according to the prescription of lAarseth et al.l (|l974l V Their 
bulk velocities are then scaled to a desired virial ratio Qi = 
r/|f2| (where T is the total kinetic energy and Q. the total 
potential energy) where Qi = 0.5 is virial equilibrium and 
our scaling ranges between Qi = and 0.5. 

Each clump is assumed to be a virialised Plummer 
sphere with a Plummer scale radius of iipi = 0.01 pc 
and a cut-off radius beyond which no stars are placed of 
-Rcut = 5-Rpi = 0.05 pc. We assume that sub-clumps are 
virialised initially as their relaxation time is so short that 
they will rapidly virialise (but this process is also effective 
at destroying clumps as we shall see). 

We take equal-mass stars of mass O.5M0 (roughly the 
average mass of a star from a standard IMF). This means 
that our clumps contain from 12 to 160 stars depending on 
the values of e and A'o. Again we acknowledge that equal- 
mass stars are not realistic. In particular, differences in stel- 
lar masses will have a significant effect on two-body encoun- 
ters which we are attempting to examine in particular de- 
tail. However, introducing a range of stellar masses would 
significantly increase stochasticity and add another free pa- 
rameter, so we choose to ignore it for now. 

Therefore the important parameters are: 

• star formation efficiency e; 

• gas Plummer scale radius _Rpi; 



• number of subclumps A*'o; 

• virial ratio of the stellar distribution Qi. 

From these parameters it is possible to calculate a number 
of very useful quantities. 

The filling factor, a, is the fraction of the volume which 
contains subclumps 



The crossing time of the whole system, T^^ , or of an 
individual clump, Tcr is the typical time taken to cross the 
whole system or individual clump (the typical size divided 
by the typical speed). 

The two-body relaxation time is a measure of how 
rapidly the internal velocities of a clump will change by or- 
der their own magnitude and is given by 

irelax = 0.1-j T-TT T^cr, (2) 
m^JVpart J 

where A'part is the total number of particles in the system, 
and tcr is the crossing-time of the system. 

A list of all the simulations and their parameters can be 
seen in Table [T] Three different random realisations of each 
parameter set are performed. 



3 RESULTS 

First we shall examine in detail our 'standard model'. This 
is a virialised Qi = 0.5 cluster with a star formation effi- 
ciency of e = 0.32 and A^o = 16. This cluster has a stellar 
mass of 32OM0 in 16 2OM0 clumps (40 equal-mass stars per 
clump). The cluster has a filling factor of a = 0.05 and a 
total crossing time of 260 kyr. Each clump has a crossing 
time of 20 kyr, and a relaxation time of 22 kyr. 

Figure [T] shows the evolution of the standard model for 
10 crossing times (2.6 Myr). Initially (t=0 Myr) the stars 
have a highly clumpy and sub-structured distribution. By 
0.5 Myr the initial structure is already less obvious, and by 
> 1 Myr the cluster has a fairly smooth appearance with 
very little evidence of the initial clumps. 

It is expected that the erasure of substructure would be 
due to one, some, or all of the following mechanisms. 

(i) Internal scattering and the ejection of stars from a 
clump by internal two-body interactions. 

(ii) Tidal stripping of clumps by the gas potential. 

(iii) Tidal encounters between clumps. 

(iv) Collisions between clumps or stars. 

An examination of Figure [T] suggests that internal scat- 
tering is a crucial factor in the erasure of substructure. 
By 0.5 Myr, the appearance of the cluster is already quite 
smooth (we shall return later to quantify the erasure of 
substructure, but for now we will use a 'by eye' examina- 
tion) . In only two system crossing times interactions between 
clumps cannot have been important. The high density of 
the initial clumps also suggests that tidal stripping by other 
clumps or the gas potential can not have been responsible 
for the smoothness. The only process that works so rapidly 
on clumps is internal relaxation as the clumps are 20 internal 
relaxation times old by 0.5 Myr. 
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t = O.OMyr 




t = 0.5Myr 




t=1.0Myi- 




t = 2.6Myr 




Figure 1. An xyz-plot of evolution of substructure in a standard model simulation. The upper- left panel shows the initial stellar 
distribution {t=0.0 Myr). Stars are initially distributed in well defined clumps. The standard model star-forming region has a crossing- 
time T^^ = 260 kyr. As the simulation evolves, we show snapshots at 2 (upper-right panel), 4 (lower-left panel), and 10 (lower-right 
panel) crossing-times. Individual clumps 'pufT-up' significantly in less than 2 crossing-times. By 4 crossing-times, only a few clearly 
defined clumps are visible, and by 10 crossing-times almost all substructure has been erased. 



The effect of internal relaxation is to eject stars from the 
clumps forming a smooth background of stars, and also to 
increase the size of the clumps ('puSing-up'). The pufBng-up 
of clumps is obvious in Figure [1] where even isolated clumps 
at 0.5 Myr are clearly larger than initially. Binary formation 
within clumps is observed. This likely plays a significant role 
in enhancing scattering and ejection of stars from clumps. 

As well as introducing a background of ejected stars, 
puffing-up has two effects which further increase the rate 
of clump mergers. As clumps are larger they are more sus- 
ceptible to tidal stripping, and the filling factor increases - 
significantly increasing the rate of clump collisions. 

In the clump merger simulations of F09 these effects 
arising from two-body encounters were missed as the SU- 
PERBOX code used for these simulations damps two-body en- 
counters entirely. This is illustrated in Figure [2] showing two 
simulations with the same initial condition^ having evolved 
for 0.7 Myr (2.7 crossing-times). Therefore it is worth revis- 

^ It should be noted that the SUPERBOX runs have 10^ much 
lower-mass particles initially in each clump and here we display 
only 40 for a fair comparison. 



iting the results of F09 when applied to star clusters in hght 
of this vitally important physical process. We note that the 
F09 results apply to situations where two-body encounters 
can be considered negligible. 



3.1 Clump counting with a Minimum Spanning 

Tree 

A minimum spanning tree (MST) is the shortest path link- 
ing a set of n points with no closed loops. The MST is useful 
in that it always has n — 1 connections (edges) wit h a un ique 
total length (see e .g. ICartwright fc WhitwoTthl (|2004 ) or 
I Allison et aj] )2009l ) for uses of the MST in a s trono my). We 
use the algorithm described in [au ison et al. I (|2009l) to con- 
struct an MST for a simulation at some point in time. 

In order to identify clumps we introduce a cutting 
length Zcut- If an edge has a length greater than Zcut then 
it is removed and we examine the subset of connections re- 
maining after cutting the longest edges. If a subset contains 
more than one-third of the stars initially within a subclump 
then we define it as a clump. (Note that iGutermuth et ahl 
(|2009l ) use a similar method, as do many friends-of-friends 
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Table 1. A complete list of the parameters of all the simulations in our parameter study. The table is split by horizontal lines into sets 
(set 1 to set 6 from top to bottom). Each set is chosen to test the influence of a specific parameter on the formation rate of the cluster 
(see text for further details). Columns give the filling factor a, the SFE e, number of clumps A^Oi followed by the Plummer radius flpp the 
cut-off radius R^'^^ , the total mass M^^ , and the crossing time of the star- forming region T^^ . The following two columns are initial Virial 
ratio Qi, and corresponding velocity dispersion of clumps with respect to their clumps within the region cr^'f^. The next two columns are 
the mass in stars A'/star and mass in gas A^fgas (modelled as an analytical background) within the star-forming region. Finally we show 
the Plummer radius _Rpi, the cut-off radius -Rcut, the mass Mpi, the crossing time Tcr, and the relaxation time of an individual clump 

^rclax ■ 

clump finders.) However, we introduce an extra element as 
our clumps are not just close in physical space, but in phase 
space. Therefore we build our MST in 6D phase space before 
applying the cut. 

A significant problem with clump finding is that it can 
be very sensitive to the cutting length used. In our initial 
conditions we have the luxury of knowing what clumps are 
present and where they are. This allows us to fine tune our 
cutting length to get the right answer initially (a bad choice 
of cutting length can result in garbage). We also check a 
number of simulations by eye to see that the structures se- 
lected as clumps are indeed clumps, and that no structure 
have been missed. This is not ideal, but the best we can do 
at the moment with no way of properly selecting a cutting 
length (|Gutermuth et al] (|2009l ) propose a way of selecting 
the cutting length, but this only works well when structures 
are distinct). 

An illustration of the method is given in Figure [3] which 
quantifies the evolution of the nbodyG and SUPERBOX sim- 
ulations illustrated in Figure [2] The solid lines shows the 
rapid decrease in the number of clumps in the nbodyG sim- 
ulation compared to those in the SUPERBOX simulation in 
the dashed line. This agrees with the by eye assessment of 



the far more rapid erasure of substructure in the former sim- 
ulation (Figure [T]). 



4 A PARAMETER STUDY OF CLUSTER 
FORMATION RATES 

With a quantitative measure of the evolution of the sub- 
structure and a better understanding and modelling of the 
physical processes behind this erasure, we can revisit the 
study of F09 to quantify the rate at which structure is 
erased. 

F09 suggested that the rate of erasure of substructure 
measured by the number of clumps at time r, N{t) com- 
pared to the initial number of subclumps A^o could be well- 
fitted by an equation of the form 

N{t) = (A^o - l)exp(-77r) + 1, (3) 

where 77 is a free parameter that depends on the initial con- 
ditions of the simulation. A large value of -q corresponds to 
a rapid loss of substructure and the rapid appearance of a 
smooth cluster. Therefore we refer to rj as the 'cluster for- 
mation rate'. 
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Figure 2. The resulting stellar distribution after 0.7 Myr of evolution from the standard model initial conditions in; an NBODY6 simulation 
(left panel), and a SUPERBOX simulation (right panel). Substructure has been erased more rapidly in the NB0DY6 simulations as a result 
of a realistic handling of two-body relaxation. 



a = 0.05,N^= 16, E = 0.32 




Crossing times elapsed x 

Figure 3. A graph of evolution of number of clumps N with 
time r (in units of the star-forming region crossing time). Both 
curves are results for simulations with the standard model initial 
conditions. The solid (black) line is the results for the NB0DY6 run 
as measured using a minimum spanning tree based clump finder. 
The dashed (red) line is the results for the SUPERBOX run. As can 
be qualitatively assessed by eye in Figure |2] substructure is more 
quickly erased in the NB0DY6 simulations as a result of a realistic 
handling of two-body relaxation. Equivalently, the NBODY6 star 
cluster forms more rapidly. 



set 4: a=0.05, e=0.10 

' I ' I ' I ' r 



simulation 

- - ■ Exponential Fit: (ri=0.39) 




Crossing times elapsed x 



Figure 4. A graph of evolution of number of clumps N with time 
T (in units of the star- forming region crossing time). The solid 
(black) curve is the results of a Set 4 simulation as measured using 
the minimum spanning tree based clump finder. The dashed (red) 
line is an example of an exponential fit, using Equation [S] to the 
simulation results. As presented in the key, this fit provides us 
with a value for the cluster formation rate parameter tj. A high 
resulting value of r] indicates that substructure has been erased 
rapidly, and thus that a smooth cluster has formed rapidly. 



As shown in Figure H] we agree with F09 tiiat an ex- 
ponential decay is indeed a good fit to the rate at which 
substructure is lost. 



4.1 Results of the parameter study on cluster 
formation rate 

In Figure [S] we show the variation of the cluster formation 
rate 77 with different parameters in our study (see Table [1]). 
Please note the differing y-axis scales in the panels of Figure 
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[5] Remember also that each clump has the same size in each 
simulation. 



enhance clump-clump collisions and tidal interactions that 
are the crucial physical parameters. 



4.2 The initial Virial ratio Q\ 

The upper-left panel of Figure [S] shows how cluster forma- 
tion rate depends on the initial virial ratio. As the initial 
virial ratio becomes increasingly sub- virial (cool), the clus- 
ter formation rate steadily increases. This occurs as a result 
of increased clump-clump interactions - for sub- virial initial 
conditions clumps tend to fall into a more compact con- 
figuration within the gas potential. Cluster formation rate 
steeply rises as Qi falls below 0.1. In this case, clumps tend 
to fall towards the centre of the gas potential well on time- 
scales ~ 1 free-fall time, resulting in multiple simultaneous 
clump-clump collisions. The initial Virial ratio Qi is a key 
parameter controlling the cluster formation rate. 



4.3 The filling factor a 

The upper-right panel of Figure [5] shows how cluster forma- 
tion rate depends on the filling factor a. For a < 0.2, the 
cluster formation rate is fairly constant as it is dominated by 
the ejection of stars from clumps rather than clump-clump 
interactions. However for very high filling factors (a ~ 0.5) 
the cluster formation rate becomes high. When filling fac- 
tor is this high, clumps are almost overlapping in the initial 
conditions, and this is enhanced by clumps pufHng-up, hence 
rapid merging occurs. 



4.4 The the initial number of clumps No and star 
formation efficiency e 

The middle-left panel of Figure [5] shows how the cluster for- 
mation rate depends on the the initial number of clumps A^o, 
and the middle-right panel shows how cluster formation rate 
depends on star formation efficiency e. At a first glance, the 
cluster formation rate appears to depend on both Ao and 
e. A trend towards increasing cluster formation rate for in- 
creasing Ao is visible in the middle-left panel. Meanwhile, 
as star formation efficiency falls, cluster formation rate in- 
creases. This trend is the reverse of what is seen in the F09 
simulations. 

However by varying Ao and e we vary the initial num- 
ber of stars in a clump. This results in a variation in the 
relaxation time treiax of any individual clump and hence the 
speed at which clumps puff-up (see Set 3 and Set 4 of Table 
[!}. As Ao increases, trciax falls. Meanwhile as e increases, 
ireiax grows. Therefore it is difficult to separate the addi- 
tional effects of varying trciax on cluster formation rate. 



4.5 A fixed clump relaxation time troiax 

In the bottom left panel of Figure[5]we show how the cluster 
formation rate varies for clumps with a fixed internal relax- 
ation time (selecting values of Ao and e such that the re- 
laxation time is 22.9 Kyr). For a constant treiax, the cluster 
formation rate is roughly constant. This strongly suggests 
that it is the internal relaxation of clumps - the ejection of 
stars to make a background and the puffing-up of clumps to 



4.6 The relaxation time of individual clumps 

This conclusion is further supported by the bottom-right 
panel of Figure [5] where the cluster formation rate can be 
seen to decrease with greater clump internal relaxation time. 
Note that these simulations are all for an initially virialised 
clump distribution, and we also exclude the a — 0.5 fill- 
ing factor simulation as the cluster formation times in these 
simulations are dominated by other effects. It appears that 
the clump relaxation time is a key parameter controlling the 
cluster formation rate. 



5 DISCUSSION 

As discussed in the introduction, star formation is a 
messy and complex process that does not initially pro- 
duce a smooth, relaxed star cluster. To give the roughly 
spherical, smooth star clusters that we often observe, 
sub-struct ure must be erased ( sec also Aarscth fc Hills 
19721 : 1 Good win 1 991; iBoilv et al.| 1999: Kroupa & Bouvier 
20031: iGoodwin fc WhitworthI |2004| : lAUison et al] l2009l : 
Fellhauer et al.ll2009l '). 

The speed at which a cluster can erase its sub- 
structure depends on its initial virial ratio (see also 
iGoodwin fc Whitworthll2004i ). but also critically on the rate 
at which sub-structure evolves internally from its initial 
state. Dense, low-A clumps have a short internal relaxation 
time and will disperse rapidly (see also iKroupa fc Bouvieil 
l2003h . Two-body encounters within a clump eject stars, 
forming a general stellar background, as well as causing 
clumps to increase significantly in size, making clump-clump 
interactions more likely, and making them more susceptible 
to tidal stripping. However, these effects are secondary to 
internal relaxation. 

By considering the middle-right panel of Figure [5] we 
can gauge the importance of clump merging in our simula- 
tions. The data points are from Set 3 and Set 4 of our pa- 
rameter study. By increasing the star formation efficiency, 
we additionally increase the mass of the clumps. For clumps 
to merge, their relative impact velocities must be of order 
the velocity dispersion of the clumps. Therefore more mas- 
sive clumps should merge more easily as seen in the cor- 
responding F09 simulations. Instead we see the opposite - 
a decreasing cluster formation rate with increasing clump 
mass. This indicates that merging plays a minor role in the 
cluster formation process. Instead it is dominated by the 
effects of inter-clump two-body encounters. For higher star 
formation efficiencies, the initial number of stars within a 
clump increases. Thus two-body encounters occur less fre- 
quently and consequently the cluster formation rate falls (as 
seen in the middle-right panel of Figure [S]) . 

The theo ry of clump merging developed in 

iFellhauer eTsU (|2002l ) and F09 is therefore not appli- 
cable in low-A (A ~ 1000) systems such as in the 
star-forming regions modelled in this study. However, in 
high-A systems such as mergers of clusters within cluster 
complexes, the effects of two-body encounters are far less 
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Figure 5. Plots of the dependency of the cluster formation rate parameter 77 (y-axis) versus the parameters investigated in the parameter 
study (x-axis); initial Virial ratio Qi, filling factor a (upper-right), the initial number of clumps A^o (middle-left panel), region star 
formation efficiency e (middle-right panel), region star formation efficiency e for fixed clump relaxation time trelax (bottom-left panel), 
and clump relaxation time treiax (bottom-right panel) . A description and discussion of these results is provided in the text (see Section 
14.11 1. Please note varying scale on y-axis. The two upper panels have a matching y-scale to each other. The four panels beneath them, 
all have matching y-scale to each other. 
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important. In these scenarios, the F09 theory should remain 
valid. 

An important conclusion to draw from this analysis is 
that clusters can change their appearance rapidly when they 
are young. A cluster that appears smooth and relaxed at an 
age of a few Myr may well not have formed that way. A 
number of authors have recently emphasised that clusters 
evolve rapidly and that current conditi ons are not always a 
good indicator of past conditions (e.g. iBastian et al.|[2008l : 
lAlhson et ai]l2009l i. 



6 SUMMARY & CONCLUSIONS 

Both observations and theory agree that stars form with 
a complex clumpy distribution within a star forming re- 
gion. Small clumps containing ~ lO's of stars form em- 
bedded within the envelope of molecular gas from which 
they formed. However, observations of star clusters of > a 
few Myr in age often show smooth, relaxed distributions 
(e.g. the Orion Nebula Cluster at ~ 3 Myr). In order for a 
smooth spherical star cluster to form, the star-forming re- 
gion must erase its initial substructure. 

We investigate the mechanisms by which substructure is 
erased by modelling star forming regions using the nbody6 
code. Our stars are initially distributed in clumps, and em- 
bedded in a static potential to mimic the gravitational influ- 
ence of the gas envelope on stellar dynamics. We conduct a 
parameter study to investigate the key parameters effecting 
the rate at which substructure is erased. Key parameters 
include the initial Virial ratio of the clumps within the gas 
potential, the filling factor of clumps within the gas, the ini- 
tial number of clumps, and the star formation efficiency of 
the total star-forming region. 

We find a number of new and different results, tha n of 
those presented in a similar study in Fcllhaucr ct al. 20091 
(F09). These differences arise predominantly due to the 
proper treatment of two-body encounters in our simulations. 
Our key results may be summarised in the following: 

(i) Clusters form predominantly from stars that are scat- 
tered out of clumpy substructure, and not by clump merging. 

(ii) As a result, the rate at which a cluster forms is a 
strong function of the relaxation time within the clumps. 
Unlike in F09, the star formation efficiency of the region 
does not effect the cluster formation rate. 

(iii) The initial virial ratio of the clumps is also a key 
parameter controlling the rate at which a cluster forms. The 
lower the initial virial ratio, the more rapidly substructure 
is erased and a cluster forms. 

As inter-clump scattering has been demonstrated to be 
of such importance to cluster formation, it is vital that mod- 
els considering the stellar dynamics of star-forming regions 
do so correctly. The use of softened gravity between stars 
will result in suppression of two-body encounters, and as 
such a key channel by which clusters form will be missed. 

Furthermore, if star clusters form by scattering of stars 
from clumps, there is an increased likelihood that a moving 
subclump can leave a trail of stars which maintain a velocity 
signature of the clump from which they originated. Similar 
trails are reported in simulations of massive stars clusters 



within a dark matter halo (|Assmann et al.ll2010l ). Such ve- 
locity structures may be observable in young clusters with 
the advent of Gaia. If so, we anticipate that these observa- 
tions could provide strong constraints on the recent forma- 
tion history of young cluste rs. We defer a det ailed study of 
this topic to a latter paper (|Smith et al.ll20l"ll in prep.). 
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